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MULTIDIMENSIONAL  INFINITELY  DIVISIBLE  VARIABLES 
AND  PROCESSES.  PART  I:  STABLE  CASE 


By 

Raoul  LePage 


1.  Introduction. 

Series  decompositions.  Involving  the  arrival  times  of  a  Poisson 
process,  have  been  given  by  Ferguson  and  Klass  [1]  for  the  non-Gaussian 
component  of  an  arbitrary  (real-valued)  Independent  Increments  random 
function  on  the  unit  Interval.  LePage  and  Woodroofe  and  Zinn  [4]  have 
rediscovered  a  variant  of  this  decomposition  in  connection  with  their 
study,  via  order  statistics,  of  the  limit  distribution  for  self-normalized 
sums  (e.g.  Students '-t) ,  when  sampling  from  a  distribution  in  the  domain 
of  attraction  of  an  arbitrary  stable  law  of  index  a  <  2. 

The  present  paper  obtains  a  characterization  of  stable  laws  on 
spaces  of  dimension  greater  than  one.  This  characterization  is  formally 
like  that  of  Ferguson-Klass  for  dimension  one,  but  with  1.1. d.  vector 
multipliers  on  the  Poisson  terms.  The  law  of  these  coefficients  may  be 
chosen  proportional  to  the  L£vy  measure,  although  this  is  not  necessary. 
These  results  take  a  particularly  elegant  form  in  the  case  of  symmetric 
stable  laws,  where  something  of  a  calculus  is  developed  showing: 

(i)  which  Levy  measure  associates  with  vector  coefficients  other  than 
the  aforementioned  ones,  (il)  what  happens  when  independent  stables  are 
linearly  combined  as  in  a  weighted  sum,  (ill)  how  to  construct  an  arbi¬ 
trary  multidimensional  Independent -Increments  symmetric  stable  set 
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function,  and  (iv)  how  to  construct  an  arbitrary  harmonizable  stationary 
symmetric  stable  random  function  having  multidimensional  domain  and/or 
range.  Symmetric  stable  laws  are  shown  to  be  mixtures  of  Gaussian  laws. 

Partly  because  of  the  self-contained  character  of  Kuelbs'  paper 
[3],  in  which  the  characterization  of  the  log-characteristic  function 
of  a  stable  law  is  extended  to  real  separable  Hilbert  space,  the 
Hilbert  space  level  of  generality  has  been  chosen  for  this  paper. 

Later  extensions  of  Kuelbs'  result  to  Banach  and  more  general  spaces 
support  a  corresponding  generalization  of  these  results.  In  addition 
to  Kuelbs'  result  we  need  a  method  employed  by  Ferguson  and  Klass  to 
transform  certain  dependent  series  into  eventually  identical  indepen¬ 
dent  ones.  We  also  require  standard  results  giving  conditions  under 
which  an  independent  series  in  Hilbert  space  converges  almost  surely 
(e.g.  [2],  Theorem  5.3).  The  rest  of  the  paper  is  basically  self- 
contained  and  affords  a  surprisingly  accessible  and  clear  view  of  a  <  2 
stable  laws,  and  random  functions,  based  on  elementary  series  constructions 

Part  II  of  this  paper  will  generalize  these  results  to  the  infinitely 
divisible  case. 


2 
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Notation 


The  following  symbols  and  conventions  will  be  in  force  throughout 
this  paper. 


(2.1) 


A 

D 


a 

{rrj  >  D 

H 


"is  asymptotic  with" 

"equals  by  definition" 

"has  the  same  distribution  as" 

"converges  in  distribution  to” 

0  <  a  <  2,  an  index  of  stability 

arrival  times  of  a  Poisson  process 
with  unit  rate 

a  real  separable  Hilbert  space 


The  material  of  the  next  section  is  drawn  from  [4], 
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3. 


Motivation. 


Limit  theorems  are  not  the  subject  of  this  paper.  However,  we 
should  not  proceed  without  benefit  of  the  following  example,  which 
exposes  some  connections  between  a  <  2  stable  r.v.  and  the  Poisson 
process. 

Let  (Cj, j  >_  1}  be  independent  of  the  sequence  T  and  i.i.d. 
with  P(e^  =  1)  =»  P(e^  =  -1)  ”  y,  and  define  G(x)  “  x-01,  Vx  _>  1  • 

Think  of  G(x)  =  P(]x|  _>  x)  Vx  1,  where  X  is  a  r.v.  symmetrically 

distributed  about  zero.  We  will  construct  particular  r.v.  X^ . X^, 

i.i.d.  as  X,  whose  normalized  sums  converge  in  distribution  to  the 
symmetric  stable  law  of  index  a.  To  do  this,  use  the  arrival  times  of  a 
Poisson  process  to  generate  uniform  order  statistics,  apply  G  1  to  these, 
multiply  by  the  signs  e,  and  permute.  As  constructed,  the  normalized 


sums  will  actually  converge  almost  surely  to  E^ 


-1/a 


(see  (3.1) 


below),  a  series  possessing  the  symmetric  stable  law  of  index  a. 

co  ~l/0t 

In  fact,  a  direct  proof  of  the  stability  of  E^  follows 

easily  from  the  observation  that  the  arrival  times  of  several  (say 
K  >  1)  independent  unit  rate  Poisson  processes  (run  simultaneously) 
constitute  K  1  times  the  arrival  times  of  a  unit  rate  Poisson  process. 
This  argument  works  just  as  well  for  e  replaced  by  any  vector  sequence 
(provided  the  series  converges)  and  suggests  the  multivariate  extensions 
of  sections  4  and  S. 

For  each  n  >  1  let  U...  <  U...  <  •••  <  U,  .  denote  the 

(1)  -  (2)  -  -  (n) 

order  statistics  of  i.i.d.  random  variables  U, ,...,U  which  are 

1  n 


uniformly  distributed  on  [0,1].  Then  for  each  fixed  n  >  1,  letting 
Xj  -  e^ay,  Vj  >  1, 


(3.1) 


-1/a  n  D  -1/a  _n  _-l,  \ 

n  ^  ^  (0(J)> 


-1/a  _n  tT-l/a 

■  ”  eju(j> 


D  -1/a  _n  .  .  -1/a 

n  Z1  Ej(rj/rn+l) 


— n  p  — l/a._  /  \  1/a 

"  21  £jrj  (rn+l/n) 


1  I 

■fa. s.  ta.s. 

oo  _-]/a  . 

Z1  1 


(SLLN) 


It  is  convenient  to  refer  to  e^,  J  >  1  as  the  residual  order  statistics, 

keeping  in  mind  that  the  ordering  is  on  decreasing  absolute  values. 

The  same  example  suggests  an  invariance  principle  (proved  in  [4]) 

for  self-normed  sums  such  as  l”  X  //z”x2  which,  regardless  of  a  , 

l  j  in 

converge  in  distribution  to  a  limit  law  depending  only  on  the  stable 
attracting  X^.  For  the  r.v.  constructed  above,  Vn  1, 


Use  Tj  ^  J  and  Z*  j“2^a  <  »  a.s.,  and  apply  the  3-series  theorem 
conditional  on  the  sequence  T. 


1. 
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(3.2) 


En  X 
L1  Xi 


ffA  F 


zne  r  -1/a 

zieJ  J... 


-2/a 


a.s. 


r“  e  j.-1/b 

a  j  3 


£ 


,-2/a 


2 

That  is,  the  limit  law  of  the  t-statistic  *  is  that  of  the  t-statistic 
calculated  on  the  residual  order  statistics  (see  also  [6]). 

Even  the  construction  of  stable  independent-increments  processes 
can  be  motivated  by  means  of  the  same  example.  We  restrict  our  atten¬ 
tion  to  the  homogeneous  increments  case.  Let  [T^  ,  j  >_  1]  2  [U^.j  >_  1], 
and  suppose  the  sequences  T,  e,  T  are  mutually  independent.  The 


partial  sum  processes  X^  ,  0  <_  t  £  1,  Vn  ^  1  can  be  effected  by 

independent  selections  of  X^,...,Xn  into  subsets  of  sizes  [nt]  using 
multiplication  by  indicators: 


(3.3) 


ljn)(t)  -  ktx  < 


ijn)(t)  -  KTj  _ 


.  ™-h<i  Ti 

n  +  1  -  j 


(n) 


(t) 


-)  ,  V  1  <  j  <  n  . 


Then  for  each  n  >  1, 


(3.4) 


(n“1/a  z[nt,XJt  te[0,l]}  2  (s"  ijn)(t)  ejrJVVn+1/n)1/o,  te[0,l]} 


+  a.s.  in  D[0,1] 


(Z*  I(Tj  <  t)ejr“1/ot,  tc [ 0, 1] }  . 

Details  of  this  argument  are  unpublished. 


2 

*  The  square  of  this  t-statistic  is  simply  related  to  the  square  of 
Students'-t,  and  both  have  the  same  limit  law. 


6 


4. 


Stable  Laws  on  H. 


'mi 


Suppose  {Xj ,j  l)  are  i.i.d.  random  vectors  in  a  real  separable 

Hilbert  space  H  and  that  the  sequences  X,  T  are  independent.  For 

each  n  ^  1  denote  by  the  number  of  arrival  times  T  in  the 

interval  [0,a  ]»  a  ■  E?  j  This  choice  of  a  is  from  [1].  Its 
n  n  l  n 

advantages  will  be  apparent  in  what  follows. 

K 

Remark.  Sums  of  the  kind  E^n  are  for  each  n  >  1  defined  to  zero  on 


the  event  K  =  0.  Use  (  ,  ), 


n 


,  to  denote  H  inner  product  and  norm. 


For  each  n  ^  1,  xeH,  c  >  0,  (see  also  [1],  pg.  1639), 


(4.1)  Ee 


K  K 

Kx.cs  nx  r:1/a)  k  icE  n(x,x ,)r~1/a 

1  i  1  «  EE  ne  1  J  J 


a  ic(x,X,)t_1/a  K 
dt)  n 


-1  fa n  „  1C'X’A1 

n  Jo 

Jo 


a  ic(x,X.)t  1^°t 
E  L  n(e  1  -l)dt 


Ecaa  J\ 


i(x,X  )r  d 
1/a  (e  1  -1) 


1+a 


Kuelbs  ([3],  lemma  2.2)  has  proved  that  the  log-characteristic 
function  of  a  (non-Gausslan)  stable  probability  measure  on  H  is 
necessarily  of  the  following  form,  for  a  unique  3  E  H  and  finite 
Borel  measure  a  on  S  ■  (xeh|||x||  -  l), 


(4.2)  i(x,g)  +  f  f  (ei(x’s)r  ~-o( ds)  ,  Vx  e  H  . 

J  S'O  l+r  r1+0t 

It  is  convenient  to  refer  to: 

(4.3)  expression  (4.2)  with  g=0. 

Define  6  ■  a(S),  )J  ■  a/6,  and 


+  +  00 
■  I L- 


sr 


1/a  fi-l/a  a-l/a  1+r2 


dr 

14a 


a(ds),  Vn  >  1 


Lemma  4.4  If  (4.3)  is  the  characteristic  function  of  a  stable  law  on 
H  and  if  the  sequence  (X^ , j  _>  l)  is  i.i.d.  y  and  independent  of  the 
sequence  T  ,  define 


x(n)  A 


a-l/<*sl/a 


y]1/a-Bn ■ 


Vn  >  1  . 


Then  VxeH,  E  exp  i(x,X  )  converges  to  (4.3). 

n 

Remark.  This  result  is  not  altogether  satisfactory  since  convergence  of 
the  series  (x^n\  n  ^  l}  is  through  stochastic  times  (K^,  n  l)  , 
and  is  not  yet  a.s.  in  H.  These  defects  are  remedied  in  Theorem  4.8 
below. 

Proof.  For  each  xeH,  n  1,  by  (4.1), 
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(n), 


(4.5)  loga  Eei(x’X  }  =  -i(x 


r*  i(x,5L)r 

^n>  +  E5  J-l/a.l/«  -!/«<•  "l> 

'a  o  a 

n 


dr 

1+a 


f  f°°  ,  i(x,s)r  .  i(x,s)r.  dr  f .  . 

I  J  -1/ct.l/a  -l/aU  1  ..2  }  1-kx  0(ds) 

•'s  ^  ot  6  a  1+r  r 

n 


(4-3)  □ 


From  ([1],  lemma  2),  we  conclude  that 


(,Kn+l  I 

j\  xirj  •  "-'r 


are  Independent.  Furthermore,  using  an  argument  drawn  from  [4],  Vn  >  1, 


(4.6)  E 


n  1 


<  a~2,a  E(K  --K  )2  ^  < 

—  n  n+1  n 


n(loge  n) 


2/a  * 


which  is  summable  in  n.  Therefore, 


Ij®  X  Tj^-EZ^  Xj(rj1/aAl) 


converges  in  probability  in  H.  A  short  calculation  gives 
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(4.7)  A  (a, 5)  5  E  a-1/o51/a  £  (r^Al)  -  .  j  _17  1/a  _1/a  -j 

J  ™  6  a  1+r  r 

n 


-a  a 


3./ 


-1/a 


,-l/afl/a  f  n  ,,,-1/a.,*  t  s 

s  Jo  <(t  dt 


finite  limit  =  A(a,6)  as  n  -*■ 


Therefore  X^n^  converges  in  probability. 

Theorem  4.8.  If  (4.3)  is  the  log  characteristic  function  of  a  stable 
law  on  H  then  the  series 


a-1  /a6 1  { X .  rT1  /<X-  (  EX. ) 

ill  1 


rj 

jj~i 


-1/a 


Ha-2/V/a  t'2/(x 


)  dt} 


converges  a.s.  in  H  to  a  random  vector  with  log  characteristic 

function  (4.3). 

Remark.  Centering  is  not  needed  for  the  case  a  <  1,  nor  is  it 
needed  for  the  symmetric  case  which  will  appear  in  [4], 

Proof.  Since  X^  converges  in  probability  and  X^  is  an  inde¬ 
pendent  series,  we  conclude  by  ([2],  Theorem  5.3(6))  that  X^n^ 
converges  a.s.  in  H.  Recall  that  with  probability  one  j  finite 
M  such  that  (n  M)  (^  smallest  N(n)  with  n  ■  K^^).  Then 
Vn  ^  m, 

~ 

The  first  term  equals  the  first  term  above.  For  the  second  term 

_  -lr  -a 
use  t  *  a  or 
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MO* 


(4 


.9)  X(N(n))  *  ct~1/0t51/c4£s^1  x  r"1/cx  -(EX  )  [  (  ^  (t~^a Al)dtl  + 

J  J  -*n 


(EX  )A  (a,6)  4 
1  n 


Since  for  n  >  M,  (n  =  ^n))=>  <  rn  1  aN(n)> 


(4.10) 


f  N(n)  _v 
I  (t  1/aAl)dt 

Jr 

n 


<  (N(n)) 


et  •  S  • 


o  . 


By  the  law  of  the  Iterated  logarithm,  a.s.  eventually  as  n  -*•  00 


(4.11) 


rrn  ,  -l/a.-v 
j  (t  A l)dt 

n 


=  f  n  t"1/0t  dt  <  |r  -n|  (n  A  r  )~1/a 
'  n 


—  n 


£  2  -/iTlog  log  n  (n  4  o(n))  0  . 


Therefore, 


(4.12)  j^(N(n) )  =  a-l/a6l/a{sn  x  r"1/a  -(EX.)  f  (t"1/a A l)dt}  - 

1  J  J  1  j  o 

A(a,6)(EX1)  4o(l)  , 


which  converges  a.s.  in  H  to  a  random  vector  with  log  characteristic 
function  (4.3).  q 

Remark.  The  centerings  used  above  also  have  an  interpretation  involving 
(EX^)  E(Fj^aAl),  which  will  not  be  given  here. 


o(l) 
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5. 


Symmetric  Case,  Multiple  Representations. 

In  this  section  we  do  not  assume  that  X^  is  distributed  according 
to  the  measure  y  ,  or  even  restrict  its  distribution  to  S.  Suppose 
X,  e,  r  are  mutually  independent  sequences,  with  e  ,  T  as  in  section 
3,  and  (X^ ,  j  1}  i.i.d.  in  H  with  e|[XjJ|  a  <  ». 

Remark.  Series  of  form  ^e  terme<*  symmetric. 

Lemma  5.1.  The  symmetric  series  converges  a.s.  in 

H  and  the  log  characteristic  function  of  its  limit  is  e] (x,X  )|aB(a), 
VxeH,  where  B(a)  =  J ^  a(cos(r)-l)  — . 

Proof.  The  arguments  needed  are  similar  to  those  of  theorem  4.8,  but 
easier  in  the  symmetric  case.  For  each  x  E  H, 

K  . 

i(x,Z  ne  X  T~  /a)(4  .  r°°  i£  |(x,X  )|r 

(5.2)  loge  Ee  1  J  3  » Ej  _1/(J  (e  -1) 

a  r 


oE  f  (e  1  -DKKx.X^I  <  raya)|(x,X1)|0t 


/  00 

aE  /  (cos(r)-l)l(!(x,X1)|  <  ray°)|(x,X1)|a 

>0  r 


•/. 


(cos(r)-l)  eKx.Xj)!01 


From  ([3],  corollary  2.1),  the  limit  in  (5.2)  is  the  characteristic 
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function  of  a  (symmetric)  stable  law  on  H.  Since  the  sums 


Z„n+1e  X  r’1^  are  independent  and  symmetric,  we  have  from  ([2], 
n  J  J  J 

K  _1/a 

theorem  5.3(1))  that  Z*  ejXjrj  converges  a.s.  in  H.  Since 

eventually  a.s.  £  K^+l  as  n  -*•  00  ,  we  conclude  ^ 

converges  a.s.  in  H.  □ 

Several  series  may  represent  the  same  stable  law. 

Theorem  5.2.  If  E||xi|j°t<  then  for  every  x  e  H, 

E|(x,X_)  |“  -  Xx |pE|  (x,X^)  1“  where  X*  is  distributed  on  S 
according  to  the  measure: 

*  x,  iix.ir 

p(xi e  A)  -  EI<pjre  A) 


EX, 


Proof.  For  every  x  e  H, 


E|(x,Xj)|a  -  e|(x,  l^ll01 

-  Ejlxjj01  Ej(x,x*)|a  .  o 

As  an  example  of  the  above,  every  symmetric  stable  law  on  H  has 
a  construction  l"  Z^  X^  rJ1^CX/E(|  Z1  {a)  1^ot  in  terms  of  an  independent 
sequence  Z  of  i.l.d.  standard  normal  r.v..  Conditional  on  the  sequences 
T  and  X  the  symmetric  stable  is  Gaussian.  That  is,  symmetric  stable 
laws  are  particular  mixtures  of  Gaussian  laws  with  zero  means  and  differing 
covariance  kernels.  The  latter  will  not  in  general  differ  only  by  scale, 
though  this  is  necessarily  true  for  H  »  R^. 
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6.  Symmetric  Stable  Random  Set  Functions. 

Let  {Tj»  j  1}  be  i.i.d.  taking  values  in  a  measurable  space 
with  measurable  sets  generically  denoted  A.  Suppose  r,  e,  X,  T 
are  mutually  independent,  where  the  latter  three  sequences  are  as  in 
section  5.  Define 

(6.1)  X(A)  -  e"  KTj  eA)X  ,  VA  . 

Theorem  6.2.  The  series  (6.1)  is  a.s.  convergent  for  each  A,  is 
jointly  symmetric  stable  for  finitely  many  A  at  a  time,  and 
X(A^)  ,...,  X(A^)  are  for  each  n  1  mutually  independent  if 
A^ . A^  are  mutually  disjoint. 

Proof .  For  each  n  ^  1,  x e  H,  real  numbers  r^,...,rn>  and  measur¬ 
able  sets  ^.....A^ 

i(x,E?  r  X(a  ))  i  E*  (E?  ,  r  I(t  eA.)}(x,X.)e.r:1/a 

(6.3)  logfi  E  e  1  k  H  „  lQg  £  g  j-1  k-1  k  j  Tc  j  j  j 

(5.1) 

-  B(a)E|(x,X1)|0EjE^1  rkI(T1eAk)|a  . 

If  A^,...,An  are  mutually  disjoint  then 

(6.4)  iICj  Vh.eVl”  '  Cl  I'kl0  P<T1E  V  '  O 
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Remark.  The  simplicity  of  this  construction  is  interesting,  as  is  the 

way  in  which  ot-dependence,  dimensional  structure,  and  functional 

dependence  are  identified  with  mutually  Independent  coefficient 
r-l/a  v  _ 

sequences  1  ,  X,  I. 

Remark.  Schilder  [ 8J  and  Kuelbs  [3]  have  explored  a  representation  of 
multidimensional  symmetric  stable  r.v.  by  means  of  a  stochastic  Integral 
with  respect  to  a  one-dimensional  stable  independent  Increments  process. 
Theorem  (4.8)  and  lemma  (5.1)  sharpen  and  extend  such  representations 
by  connecting  them  with  the  Ferguson-Klass  representation,  making 
explicit  the  choice  of  coefficients  required  to  obtain  each  stable 

4. 

law  *,  and  establishing  H  convergence  of  the  indicated  series. 

oo  —1/a 

Remark.  Suppose  K  >  1  and  -  Z^  ^kj^kj  *  *  —  ^  —  K»  are 
independently  constructed  (as  per  (5.1))  symmetric  stable  r.v.  taking 
values  in  H.  Then  for  an  arbitrary  choice  of  real  numbers  r^,...,rR 
the  sum  E^  r^Yfc  is  representable  E^  where  j  _>  1} 

are  1.1. d.  and  -  an  equiprobable  random  selection  from  *  Vu- 

1  /tt 

...,K  '  rn^n2*  This  uses  the  property  (discussed  in  section  3)  of  K 
Poisson  processes  run  simultaneously. 


4 

Including  the  non-syametrlc  stable  laws. 
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7.  Harmonizable  Stationary  Symmetric  Stable  Random  Functions. 


Basically,  we  seek  to  construct  the  stable  analogues  of  Gaussian 
stationary  random  functions  having  a  harmonic  decomposition.  The 
characteristic  function  of  such  a  Gaussian  random  function  involves 


c,  ‘k  ‘l(Vtk>° 


where  t  is  generic  for  a  point  of  the  domain,  and  (A  ,  )  is 

1  o 

a  random  linear  function  on  the  domain.  The  stable  analogues  of 

these  Gaussian  random  functions  have  characteristic  functions  that 
th 

employ  an  a  -power  in  this  integral  instead  of  the  2,  but  are 
otherwise  the  same.  Define  Vt, 


(7.1) 


X(t)  £  z"  cos((Aj>t)o+  ®j>xjejrj1/o 


where  A  are  i.i.d.,  0  are  i.i.d.  uniforms  on  [-11,11],  {E,X,r) 
are  as  in  section  6,  and  A,  0,  c,  X,  T  are  mutually  independent 
sequences.  The  series  (7.1)  is  a.s.  convergent  in  H  for  each  t 
by  lemma  S.l.  The  random  function  X(  )  is  clearly  stationary 
because  6  are  uniform  on  [-11,11],  but  this  will  also  be  a  simple 
consequence  of  the  form  of  the  characteristic  function  which  we  now 
compute. 


(7.2) 


lo«e 


E  e 


Kx.Zk-i  rkX(tk^ 


1  OCl  r 

-  log#  Ee  3  1 
(5i1)  B(a)E|(x,X1)|a  E|E 


k  cos((AJ,tk)o40j)}ej(x,Xj)r‘1/° 
W  rk  co‘«'Wo  +0!>>a  • 
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The  final  term  in  the  right  side  of  (7.2)  reduces  as  follows,  with 


z  A  n  ^Wo 
Z  E1  rk  6 


i0.  -10 


(7.3)  E|Z“  rk  cos((A1,tk)o  +01)|“  -  E2^|ze  1  +  ie 


l|Ct 


-210 

-  e|z  |0t2-0tEz  1 1  +  z-1ze  1|Ct 


«  e| z |a  2"aC(a) 


where 

c(a) .  r  ii+einia*i-  r 

j-V  J- TT 

for  all  real  i|>.  We  have  therefore  proved. 


Theorem  7.4.  The  random  function  defined  by  (7.1)  converges  a.s.  in 
H  for  each  t,  and  has  log  characteristic  function 


2'<lB(a)C(a)E|z^1 


rke 


i(/Wo 


a 


for  all  n  >  1,  r, , . . . ,r  ,  t, , . . .  ,t  .  n 

Corollary  7.5.  The  random  function  (7.1)  is  non-ergodlc  for  each  a  <  2. 
Proof.  By  using  (7.2)  the  construction  (7.1)  rmsains  valid  if  Ej  are 
replaced  by  Zj/(E|ZjJa)^a,  j  >_  1,  where  Z  is  an  Independent  i.i.d. 
standard  normal  sequence.  Conditional  on  the  sequences  A,  X,  T,  the 


process  (7.1)  Is  a.s.  stationary  Gaussian  with  discrete  spectrum 
therefore  conditionally  non-ergodic  a.s.. 


8.  Operator  Stable  and  Infinitely  Divisible  Multidimensional  Laws. 

Infinitely  divisible  laws,  of  which  the  operator  stable  laws  are  a 
special  case  with  particularly  interesting  structure,  are  treated  in 
Part  II.  In  brief,  this  is  what  happens:  A  construction  of  infinitely 
divisible  random  vectors  is  given  by  E{XjH(Tj ,X^ )-y  } ,  in  which  the  real 
function  H  is  monotone  decreasing  and  positive  for  each  value  X^ ,  and  is  deter¬ 
mined  from  the  Levy  measure.  A  construction  of  full  operator  stable  random 
vectors  in  a  finite  dimensional  real  vector  space  Rd  is  E{A(r^)Xj-Yj }»  in 
which  the  vectors  y  are  non-stochastic  centerings,  (A(t)  *  exp(Blog  t),t  >  0} 
is  the  group  of  linear  transformations  figuring  in  the  definition  of  operator 
stability  (e.g.  Sharpe  l9|),  and  the  vectors  X  are  i.i.d.  from  a  probability 


measure  (a  factor  of  the  Levy  measure)  on  a  set  of  generators  of  the  subgroups 
induced  by  A.  The  methods  of  sections  4  and  5  carry  over,  as  will  now  be  in¬ 
dicated.  If  X  is  any  i.i.d.  sequence  in  Rd,  and  X  is  independent  of  T,  then 
VxfRd,  n  >  1, 

E(ei(x,A(t)X1)_  d|  _ 

_n  t 

As  usual,  the  symmetric  case  is  simplest.  If  we  examine  Sharpe's  Theorem  5, 
we  discover  that  the  limit  of  (8.1)  is  precisely  the  form  taken  by  the  oper¬ 
ator  stable  in  this  case,  provided  we  choose  for  the  distribution  of  X^  the 
probability  measure  figuring  in  Sharpe's  representation  of  the  Levy  measure 
as  a  mixture,  this  measure  being  placed  on  (Sharpe's  notation)  generators 

0  characterized  by  s  Mg(t  9  :t>s}_s,Vs>0.  Arguing  as  in  section 
oo-l  d 

5,  we  conclude  E^  A(Tj  )X^  converges  a.s.  in  R  and  has  the  log-characteris¬ 
tic  function  which  is  the  limit  of  (8.1). 


(8.1) 


log  E  e1^'1!"  A(ri  )XJ> 


A 

^  a 
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9.  Priority  of  P.  Levy. 

P.  Levy  has  anticipated  the  series  constructions  of  one  dimensional 

stable  r.v.  with  a  <  2.  For  the  case  of  a  positive  stable  with  a  <  1, 

00 

up  to  scale  and  location,  this  construction  is  J  T.  ,  with 

j=1  J  + 

(1^,  j  >  1}  being  the  arrival  times  of  a  Poisson  process  (on  R  )  having 
unit  intensity  function.  Levy  writes  the  series  in  the  form  J  U  ,  where 


(9.1) 


(Ux>  x  >  0}  are  independent  r.v.  and 


(9.2)  P(Ux  =  x)  *  l-P(Ux  =  0)  . 

x 

Here  is  my  abstract  of  the  key  parts  of  Levy's  (1935)  arguments  for  the 
above  case: 


(9.3)  [ 


adx 


-a  ^ 

1-kx  *0 


[{U 


Ux  *  0,  x  > 


xo} 


is  finite  for  x^  > 


0] 


and  also. 


(9.4) 


■*0  adx 
0  xa 


a  „(l-a) 

1-a  0 


+  0]  * 


tE  y  u 

/.  x 

x<  x„ 


+  0] (as 


+  0)  . 


Therefore,  for  arbitrary  c^,  c^  >  0  (defining  c^  *  c®  +  Cj  and  taking 
independent  copies) , 
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(9.5)  c.  I  U<»  +  c,  l  U<2>  7  Y(1>  +  Jy<2> 

1  ^  x  2  L  x  '•  x  L  x 

XXX 


Dist. 


It 


(x) 


Dist . 


c^  \  (*>  stable) , 

x 


/ cl  i 

where  (Yv  ' ,  x  >  0}  have  respective  intensities  c,  a  dx/x  ,  k*  1,2,3, 

X  K. 

and  are  Independent  for  k  ■  1,2. 

The  above  arguments  do  yield  a  proof  of  the  representation  if  we 

apply  them  to  the  independent  sub-sums  xe  [bn+^,  n  — 

-a 

where  b^  *  log  n.  This  is  essentially  the  argument  of  Ferguson-Klass 

(1972) .  The  particular  choice  of  bn>  n  1  is  one  which  ensures  that 

eventually  as  n  <»  each  sub-sum  contains  at  most  one  summand,  so  it 

really  is  (almost)  as  though  one  could  add  independent  Ux  one  at  a  time 

toward  x  -*■  0.  A  quite  different  justification  is  to  interpret  £  Ux 

x 

as  a  generalized  process  driven  by  "white  noise"  (Ux,  x  >  0} . 

Levy's  observations  are  easily  overlooked.  Ferguson-Klass,  Vervaat 
(1979),  LePage-Wood roof e-Zinn  (1979)  (in  manuscript  form),  rediscover 
the  Levy  construction  as  byproducts  of  the  following  independent  pursuits 
respectively:  (F-K)  -representing  the  positive  non-Gaussian  part  of  an 
independent  Increments  random  function  as  the  sum  of  its  ordered  jumps. 
(V)-examining  a  shot-noise  associated  with  the  asymptotic  behavior  of  the 
solution  of  a  stochastic  difference  equation  as  time  is  increased, 
(L-W-Z)-«tudying  the  limit  behavior  of  the  normalized  order  statistics 
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from  a  distribution  attracted  to  a  stable.  Resnick  (1976)  reconciles  the 


Ferguson-Klass  construction  with  the  Ito  representation,  meaning  by  the 
latter  Ito's  generalization  of  Levy's  stochastic  integral  construction 
by  a  Poisson  random  measure. 
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